Plasmons in disordered nanoparticle chains: Localization and Transport 
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Disorder-induced effects on plasmon coupling in chains of metallic nanoparticles are studied within 
a dipole model, by considering two types of disorder: fluctuations of the particles' shapes and fluc- 
tuations of their positions. Typical localization effects are found both in the eigenmodes and in the 
transport behavior of the system, and an estimate of the localization length is made. It is argued 
that chains with deliberately introduced disorder constitute promising systems for studying localiza- 
tion effects of electromagnetic waves at optical frequencies under well controllable and manipulable 
conditions. 
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I. INTRODUCTION 

Disorder-induced localization was studied in 
detail for the wavefunctions of electrons by 
P. W. Anderson^. This localization effect rests 
on the interference of multiply scattered waves. 
Therefore, a similar effect is also present for 
electromagnetic waves^*^, with the difference 
that there is no interaction between the pho- 
tons and, hence, localization does not have to 
compete against interaction effects. In gen- 
eral, in three dimensions Anderson-localization 
of waves is observable if the density of scat- 
terers exceeds a critical value, as expressed by 
the loffc-Rcgel criterion^: kl < 1, with k de- 
noting the wavenumber, and I the mean free- 
path length. Hence, for electromagnetic waves 
in the visible spectrum (with k on the order of 
10^ m~^) the mean free-path length has to be 
on the order of a few nanometers. This can 
be achieved, for example, with semi-conductor 
powders^. 

In the present paper, localization of an elec- 
tromagnetic excitation in disordered chains of 
metallic nanoparticles is studied. In the context 
of integrated optics and plasmonic waveguides, 
the response of nanoparticle chains without 
sizeable disorder to an optical excitation with 
a frequency close to the plasmon resonance fre- 
quency of a single chain member has already 
been investigated in detail both experimen- 
tally^i^ and theoretically2r2£. In these studies 
it has mostly been assumed that the nanoparti- 
cles are perfect spheres, and it has been shown 
that such ideal chains offer interesting perspec- 
tives for sub- wavelength waveguiding. 

Recently, Alii and Engheta have made an im- 
portant step by studying the effect of small ran- 
dom uncontrollable disorder on the guidance 
properties of such chains, within first order per- 
turbation theoryi^. These authors were able 
to quantify the resulting additional radiation 
losses for the guided mode. Here, a further step 



is taken by evaluating even strong disorder, and 
making the connection to localization physics. 

In a wider context, localization of surface 
plasmons has already been reported experi- 
mentally for semicontinuous metal filmsi^, and 
studied theoretically, e. g., for randomly dis- 
tributed metal particles^. Also the celebrated 
surface enhanced Raman spectroscopy (SERS) 
is basically enabled by the formation of "hot 
spots" near a rough surfacei^. These hot spots 
are regions with an extremely enhanced near- 
field, and result from the localization of surface 
plasmons. 

The disordered nanoparticle chains studied 
in this contribution offer an alternative ap- 
proach to such phenomena. When the individ- 
ual nanoparticles are subjected to shape fluc- 
tuations, and randomly oriented, the polariza- 
tion tensors of the chain elements have, in gen- 
eral, three distinct principal axes, which are 
not aligned within the chain. As will been 
shown by studying their optical properties, such 
chains exhibit typical signatures of localized 
modes, when excited at frequencies close to the 
plasmon-resonance frequency of a single parti- 
cle. In contrast to systems with truly random 
disorder, such systems may allow one to intro- 
duce specifically engineered types of disorder, 
and therefore enable one to create hot spots at 
desired positions, and with predefined proper- 
ties. 

The model used to describe the optical prop- 
erties of the chains is briefly introduced in 
Sec. ini In Sec. IIIII numerical results for two 
types of disorder are shown and some conclud- 
ing remarks are given in Sec. IIVI 



II. MODEL 

The present theoretical study is based on a 
model of electromagnetic transport through a 
chain of metallic nanoparticles which regards 
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the chain as a hnear array of point dipolesSr^l. 
This approximation is viable as long as the in- 
dividual chain members are small compared to 
the relevant wavelengths, so that their optical 
response can be described in the quasi-static 
approximation. Moreover, it is required that 
the distance between the chain elements is large 
compared to their linear size, so that higher 
multipoles are negligiblei^. Experiments con- 
ducted in this parameter regime are in good 
qualitative agreement with predictions made by 
this model^. This dipole model was proposed in 
2000 by Brongersma et ali^, and since then has 
been elaborated further in several studies. In 
particular, Weber and Fordid pointed out that 
while in the parameter regime used, the quasi- 
static approximation is valid for the response of 
a single chain member, for the coupling of the 
particles the quasi-static model is not sufficient. 
Therefore, in the following always the standard 
field produced by an oscillating dipole pe'"* is 
used^a 
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with the vector r pointing from the source 
to the observation point (r denotes the cor- 
responding unit vector) and the wave-number 
k = uj/c. The time-dependence is suppressed 
in Eq. ((TJ and in the following (always assum- 
ing a harmonic time-dependence e"^*). 

Assuming, to begin with, that all particles 
are spheres surrounded by vacuum, the dipole 
moment p{rj) induced in the j-th particle lo- 
cated at the position rj is related to the electric 
field E{rj) through its effective polarizability 

OieS- 
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with the vacuum permittivity eq. The effective 
polarizability, in its turn, is related to the quasi- 



static one a vie 
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with k — bj / c and c denoting the vacuum speed 
of light. Here, the radiation damping is in- 
cluded within the Abraham-Lorentz approxi- 
mation^^. Therefore, the model can be ex- 
pected to be quantitatively correct as long as 



radiation damping is a small correction. Since 
the electric field at each chain site is the sum 
of the applied external field and the fields pro- 
duced by all other dipolcs, Eq. © describes a 
coupled system of equations for the moments at 
the sites. 

When replacing the spheres by arbitrarily de- 
formed particles their polarizabilities have to be 
described by tensors, instead of scalar polariz- 
abilities a and aeff. Then the induced dipoles 
and the local electric fields are no longer par- 
allel. Accordingly, the effective polarizabilities 
then are defined by diagonalizing that tensor, 
and using Eq. ([3]) for each diagonal element. 
The tensorial character of the polarizabilities 
implies that a decoupling of transversal and lon- 
gitudinal modes is no longer possible in chains 
of randomly oriented, distorted particles. If a 
chain of ideal spheres is excited with a field po- 
larized perpendicular to the chain-axis, symme- 
try requires that also all induced dipoles are 
oriented perpendicular to this axis (transver- 
sal mode); similarly a parallel excitation leads 
to dipoles oriented along the chain axis (lon- 
gitudinal mode). While this categorization is 
very helpful for modeling the perfect chain and 
also enables one to derive some analytical re- 
sults in this case, it is no longer valid for a dis- 
ordered chain. Of course, the discrimination of 
transversal and longitudinal excitation remains 
meaningful. 

This set of equations ©, augmented by ten- 
sorial polarizabilities for each particle, allows 
one to investigate the propagation of an op- 
tical excitation through a chain of disordered 
metallic nanoparticles in detail. In the present 
study mainly chains of randomly oriented el- 
lipsoids are considered, and both eigenmodes 
and transport properties are calculated. As a 
second type of disorder also chains of perfect 
spheres but with variations of the particle po- 
sitions is used. 



III. NUMERICAL RESULTS 

Let us start by investigating a chain of 50 
ideal spheres, with a center-to-ccnter distance 
of 80 nm and a sphere-radius of ao ~ 25 nm, 
assuming that the dielectric properties of the 
particles are described by the Drude model 



e(w) = eo 



with plasma frequency Wp = 1.4 ■ lO^^s"^ and 
relaxation time r = 1.27 • 10~^^s, as appropri- 
ate for silver—. In Fig. [T] an eigenmode^^ of this 
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chain is depicted by plotting the induced dipole 
moments \pn\ = |p(r„)|; this mode can be rec- 
ognized as a standing wavoii. This changes 
drastically when one individual sphere in the 
chain center is replaced by an ellipsoid, here 
with semi-axes 0.98 ao, 0.96 ag, and 1.02 aq, rep- 
resenting an isolated defect. As shown in Fig.[T] 
this defect leads to the emergence of an expo- 
nentially localized mode, similar to the localiza- 
tion of electronic wavefunctions around isolated 
lattice defects. 
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Figure 1. (Color online) Extended, standing-wave 
like mode in a chain of 50 silver nanoparticles con- 
sisting of perfect spheres (dashed line with circles) 
compared to a localized eigenmode which emerges 
if one single sphere is replaced by an ellipsoid (solid 
line with crosses). The resonance frequency of 
this mode is oj = (1.0912 — 0.147296i) ojspp, where 
ojspp = u)p/^/i is the resonance frequency of a 
sphere. 

Here, it should be mentioned that Anderson 
localization is not the appearance of a local- 
ized eigenmode in a system with one isolated 
defect, but, as already indicated by the title of 
the original worki, the absence of diffusion in 
a system containing many defects. The collab- 
orative effect of these defects leads to a break- 
down of the diffusive transport and hence, a 
change in the transport characteristic from an 
algebraic behavior, typical for a diffusive trans- 
port, to an exponential one is a clear indica- 
tor for Anderson localization. Certainly, in the 
context of electromagnetic waves also absorp- 
tion could lead to an exponential transport be- 
havior and therefore, to make the connection 
between Anderson localization and the propa- 
gation of electromagnetic waves through a dis- 
ordered medium it is important to ensure that 
the exponential decay is caused by the disorder 
and not by absorption processes. 

For studying Anderson localization in chains 
of nanoparticle, the transport through 80- 
element chains of randomly oriented ellipsoids 
with normally distributed semi- axes is investi- 
gated, assuming a mean of ao and a standard 



deviation of aa^. In Fig. [5] the results obtained 
when only the first particle is subjected to the 
incident field is plotted for a = 0.05, and com- 
pared with the transport through a chain of 
spheres. Again the center-to-center distance 
was 80 nm and oq = 25 nm; here and in the 
following the results in the disordered cases are 
averaged over 1000 realizations. 
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Figure 2. (Color onhne) Induced dipole moments 
in chains of 80 randomly oriented ellipsoids and 
ideal spheres, respectively, when the incident field 
is present only at the position of the first particle, 
possessing the frequency ui = uOspp- In a) the ex- 
ternal field is polarized perpendicular to the chain 
axis, whereas it is polarized parallel to this axis 
in b). The dashed line in a) shows that here the 
dipole moments decay almost exponentially in the 
transversely excited disordered chain. 

In the case of transversal excitation shown in 
Fig. [5] a) the disorder obviously causes a sub- 
stantial increase of the induced dipole moments 
as compared to the ideal chain; observe the log- 
arithmic scale. Most significantly, the decay of 
the moments along the chain changes from alge- 
braic to nearly exponential when disorder is in- 
troduced, as indicated by the dashed hue. This 
is, as explained above, a genuine hallmark of 
disorder-induced localization. In contrast, for 
longitudinal excitation the disorder only leads 
to a marginal decrease of the induced dipole 
moments. 

Crucial for the considerable effect of the dis- 
order in the case of transversal excitation is that 
the disorder disturbed the symmetry of the sys- 
tem. Hence, the effect does not depend on the 
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precise type of disorder, only the parallelism of 
the incident field and the induced moments has 
to be circumvented. The strong effect induced 
by this rotating of the dipoles results from the 
heavy direction-dependency of the dipole field 
(see Eq. ©). 
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Figure 3. (Color online) Induced dipole moments in 
chains of 80 spheres. In the ideal case the particles 
are arranged on a ID linear chain with a center- 
to-center distance of 80 nm, while the disorder is 
introduced due to normally distributed variations 
of the particles' positions (with a mean of and 
a standard deviation of 5nm). The positions are 
varied in all directions [disordered, see Eq. ©), or 
only in one dimension (perpendicular (see Eq. ([6)l) 
or parallel (see Eq. (0) to the incident field), re- 
spectively. The incident field is present again only 
at the position of the first particle, possessing the 
frequency uj = ujspp and polarized perpendicular 
to the chain axis. For clarity in the ID cases only 
the induced moments in every third particle are de- 
picted. 

This is demonstrated in Fig. |3] where the dis- 
order is not introduced due to distortions of the 
elements' shapes, but due to variations of the 
positions of the particles. If one assumes that 
the chain axis is oriented along the x-axis of 
the coordinate system used the sites of the ideal 
chain could be labeled by 

rjf = Tide, (4) 

with the ccnter-to-ccntcr distance d. The dis- 
order is now introduced due to normally dis- 
tributed deviations ajj'^'^ from these ideal posi- 
tions 

r*'' = nfie, + <;e, + <jey + a^e;,. (5) 

In Fig. [3] the results for three cases of this type 
of disorder are depicted for a transversal exci- 
tation. In every case the mean of a^'^'^ was 
zero with a standard deviation of 5nm and a 
site distance d of 80 nm. This disorder again 
leads to an exponential decay as demonstrated 
by the curve labeled disordered in Fig. [31 How- 
ever, for this type of disorder it is easily possible 



to change the disorder such, that the induced 
dipole moments and the incident field remain 
parallel. Due to the structure of the dipole field 
(see Eq. ([1])) a variation of the chain elements' 
position perpendicular to the incident field does 
not affect the orientation of the induced dipoles. 
To clarify this, assume a polarization of the in- 
cident field in z-direction with the chain-axis 
again in the a;-direction. Since, for a perfect 
sphere, the induced dipole moment and the field 
at the position of the sphere are parallel, all in- 
duced moments in the chain and the incident 
field are parallel if the fields at all chain sites 
are parallel to the external field. If now the 
positions of the particles are varied only per- 
pendicular to the polarization of the external 
field, for example only in x-dircction, 

r"^^' ^ jide^ + ale^, (6) 

for all 11, the dot product of r*'' with the polar- 
ization of the incident field (oc e^) vanishes, and 
the component of the field along r in Eq. ([1]) 
is zero. Therefore, the fields produced by the 
dipole moment in z-direction in the first par- 
ticle, caused by the external field, are also po- 
larized in z-direction at all positions rf^** (see 
Eq. ([IJ). Since, similar arguments hold for the 
fields produced by the dipole moments of all 
other chain members, all induced dipole mo- 
ment arc parallel to the incident field. Thus, 
a disorder described by Eq. ^ changes the 
induced moments only marginally, as demon- 
strated by the red crosses in Fig. |31 The same 
qualitative result is obtained if instead of non- 
zero a^, non-zero aj^ or even non-zero and 
non-zero are used. But, if the positions are 
varied in z-direction 

r'^i' = nde^ + ale^, (7) 

the dot products of r with the polarization of 
the incident field vanish no longer and hence, in 
this case the induced dipole moments and the 
external field are in general not parallel. This 
results in a qualitative change of the transport 
behavior as exemplified by the gray circles in 
Fig. El 

Naturally, the transport characteristics of the 
disordered chain depend on the strength of the 
particular kind of disorder, in addition to the 
center-to-center distances and to the optical 
properties of the nanoparticles. For studying 
the dependence of the distortion strength again 
chains of randomly oriented ellipsoids with nor- 
mally distributed semi-axes with a mean of 
flo and a standard deviation of aa^ are used. 
In Fig. m the transport is depicted for dif- 
ferent strengths aa^, for the interesting case 
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of transversal excitation. Here, the cross-over 
from an algebraic to the exponential decay is 
visible. In order to extract the dependence 
of the localization length b on the distortion 
strength a, the dipole moments in the middle 
part of the chain have been fitted to the expres- 
sion 



\Pn 



a exp(— 



(8) 



The results are shown in Fig.|l}D). As expected, 
an increase of the disorder strength causes a de- 
crease of the localization length, and hence a 
"stronger" localization. Since Anderson local- 
ization is caused by multiple scattered waves, 
the variation of the individual chain members 
should not be too large. After increasing the 
value of a to 0.2, for instance, one does no 
longer find an exponential decay of the induced 
dipole moments along the chain. For the ran- 
domly oriented ellipsoids this is due to the fact 
that the chain is excited with the resonance 
frequency of an undistorted chain member and 
that the distortion of the particles also affects 
the resonance frequencies; a too large distor- 
tion decreases significantly the scattering cross- 
section of the particle at the frequency used and 
thus decreases the effect of multiple scattered 
waves. 

Perfect exponential localization would imply 
that for longer chains the induced dipole mo- 
ments can drop substantially below those in 
ideal chains. Somewhat surprisingly, this does 
not seem to occur in the system studied here. 
Instead, in longer chains there appears to be 
a second cross-over after which the response of 
the disordered chain again approaches that of 
the ideal one. This is exemplified in Fig. [5] 
for a chain of 200 particles. Thus, there are 
at least two characteristic length scales: That 
which marks the onset of the fully developed 
exponential decay, and that which marks the 
return to the behavior of the ideal chain. 



particularly high field enhancement for applica- 
tion purposes. 

a) 1 




Figure 4. (Color online) Dependence of the trans- 
port characteristics on the strength a of the distor- 
tion of the splieres, for transversal excitation. Panel 
a) shows the induced dipole moments for different 
values of a, while panel b) depicts the dependence 
of the localization length b (in numbers of particles, 
see Eq. ((8}) on cr, with numerical errors as indicated 
by the symbols. In all cases the frequency of exci- 
tation was UJ = UJavv 
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IV. CONCLUSIONS 

In conclusion, chains of randomly or even 
deliberately distorted nanoparticles constitute 
novel experimentally accessible and theoreti- 
cally tractable systems for studying localiza- 
tion effects of electromagnetic waves at optical 
frequencies. Utilizing these localization effects 
in combination with specific types of the engi- 
neered disorder offers the prospect of obtaining 



Figure 5. (Color online) As Fig. [5^), but for a chain 
of 200 particles. Note that the induced dipole mo- 
ments in the disordered chain do not drop below 
those of the ideal one. 
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